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ABSTRACT 

We study the observational signatures of a potential population of low-luminosity 
quasars at high redshifts in a ACDM cosmology. We derive the evolution of the 
quasar luminosity function at fainter luminosities and higher redshifts than currently 
detected, based on three assumptions: (1) the formation of dark-matter halos follows 
the Press-Schechter theory, (2) the ratio of central black hole mass to halo mass is the 
same for all halos, and (3) the light-curve of quasars, in Eddington units, is universal. 
We show that a universal light -curve provides an excellent fit to the observed quasar 
luminosity function at redshifts 2.6 < z < 4.5. By extrapolating the evolution of this 
luminosity function to higher redshifts (4.5 < z < 20), we find that the associated early 
population of low-luminosity quasars reionizes the universe at a redshift z ~ 12. The 
reprocessing of the UV light of these quasars by dust from early type II supernovae, 
distorts the microwave background spectrum by a Compton y-parameter, y ~ 10~ 5 , 
comparable to the upper limit set by COBE. The Next Generation Space Telescope 
could detect tens of quasars from redshifts z > 10 per square arcminute, with its 
proposed InJy sensitivity at 1-3.5/mi. Absorption spectra of several such quasars 
would reveal the reionization history of the universe. 
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1. Introduction 

One of the outstanding problems in cosmology is the nature of the first generation of 
astrophysical objects which appeared when the universe first transformed from its initial smooth 
state to its current clumpy state. Although we have observational data on bright quasars and 
galaxies out to redshifts z ~ 5 (Schneider, Schmidt & Gunn 1991; Franx et al. 1997), and on 
the linear density fluctuations at redshift z ~ 10 3 (Bennett et al. 1996), there is currently no 
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direct evidence as to when and how the first structures formed, and what kind of objects were 
responsible for end of the "dark age" of the universe (Rees 1996). 

Popular Cold Dark Matter (CDM) models for structure formation predict the appearance of 
the first baryonic objects with masses M ~ 1O 5 M0 at redshifts as high as z ~ 30; objects with 
progressively higher masses assemble later (cf. Haiman, Thoul & Loeb, and references therein). 
Following virialization, the gas in these objects can only continue to collapse and fragment if it 
can cool on a timescale shorter than the Hubble time. In the metal-poor primordial gas, the only 
coolants that satisfy this requirement are neutral atomic hydrogen (H) and molecular hydrogen 
(H2). However, H2 molecules are fragile, and are easily photodissociated throughout the universe 
by trace amounts of starlight (Haiman, Rees, & Loeb 1996, 1997; Gnedin & Ostriker 1997). Hence, 
most of the first stars are expected to form inside objects with virial temperatures T v ; r ^ 10 K. 
Depending on the details of their cooling and angular momentum transport, the gas in these 
objects is expected to either fragment into stars, or form a central black hole exhibiting quasar 
activity. Conversion of a small fraction (~ 1 — 10%) of the gas into stars or quasars could reionize 
the universe, and strongly affect the entropy of the intergalactic medium. 

In two previous papers (Haiman & Loeb 1997a, hereafter HL97; and Loeb & Haiman, 1997, 
hereafter LH97), we explored the impact of the first stars on the reionization history of the 
universe. In this paper, we study the related signatures of early quasars, and compare their effects 
to those expected from the stars. For both types of sources, we calibrate the total amount of 
light they produce based on data from redshifts z <; 5. The efficiency of early star formation 
is calibrated based on the observed metallicity of the intergalactic medium (Songaila & Cowie 
1996, Tytler et al. 1995), while the early quasars are constrained to match the quasar luminosity 
function at redshifts z <; 5 (Pei 1995). In HL97 and LH97 we have studied a range of standard 
Cold Dark Matter (CDM) cosmologies, but in this paper we focus on a particular cosmological 
model with a cosmological constant, namely the "concordance model" of Ostriker & Steinhardt 
(1995). Within this model, we predict the number of high-redshift low-luminosity quasars using 
the Press-Schechter formalism (Press & Schechter 1974). 

We assume that the luminosity history of each black hole depends only on its mass, and 
postulate the existence of a universal quasar light-curve in Eddington units. This approach 
is motivated by the fact that for a sufficiently high fueling rate, quasars are likely to shine at 
their maximum possible luminosity, which is some constant fraction of the Eddington limit, for 
a time which is dictated by their final mass and radiative efficiency. If the radiative efficiency 
were determined by a universal accretion geometry, then the bright phase of quasars would 
admit such a universal light-curve. We also assume that the final black hole mass is a fixed 
fraction of the total halo mass, and allow this fraction to be a free parameter. Based on this 
minimal set of assumptions, we demonstrate that there exists a universal light curve, for which 
the Press-Schechter theory provides an excellent fit to the observed evolution of the luminosity 
function of bright quasars between redshifts 2.6 < z < 4.5. Furthermore, our fitting procedure 
results in black hole to halo mass ratios close to the typical values found in local galaxies 
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(Kormendy et al. 1997; Magorrian et al. 1997). Given this ratio and the fitted light-curve, we 
then extrapolate the observed LF to higher redshifts and low luminosities. 

Previous work on modeling the evolution of the quasar LF has been mostly phenomenological, 
and involved fitting ad-hoc parametric functions to the observations (Pei 1995, Boyle et al. 1991). 
Small & Blandford (1992) had used a universal quasar light-curve, together with the observed 
LF evolution, to infer the birth-rate of quasar black holes. Haehnelt & Rees (1993; see also 
Efstathiou & Rees 1988) adopted the Press-Schechter theory for halo formation, specified the 
quasar light-curve a priori, and introduced a redshift and halo-mass dependence to the black 
hole formation efficiency, in order to fit the observed evolution of the quasar LF. The postulated 
(ad-hoc) exponential decrease in the efficiency for black hole formation in small halos (with 
circular velocities less than 400 km/s) could be attributed to the expulsion of gas from shallow 
potential wells by supernova-driven winds or mergers. However, this prescription appears to be 
violated in nearby low-mass galaxies, such as the compact ellipticals M32 and NGC 4486B. In 
particular, van der Marel et al. (1997) infer a black hole mass of ~ 3.4 x 10 6 M Q in M32, which is 
a fraction ~8x 10~ 3 of the stellar mass of the galaxy, ~ 4 x 1O 8 M , for a central mass-to-light 
ratio of 7y = 2; while Kormendy et al. (1997) infer a black hole mass of 6 x 10 8 M Q in NGC 
4486B which is a fraction ~ 9% of the stellar mass. The efficiency of black hole formation in 
these systems had been comparable if not higher than in galaxies with deeper potential wells. 
Observations of galactic nuclei in the local universe imply black hole masses which are typically a 
fraction ~6x 10~ 3 of the total baryonic mass of their host (see summary of dynamical estimates 
in Fig. 5b of Magorrian et al. 1997; and the maximum ratio between a quasar luminosity and 
its host mass found by McLeod 1997). We therefore adopt the simplest possible assumption, 
namely that the black hole to halo mass ratio is constant, and infer the quasar light -curve from 
observations rather than postulate it. The actual formation process of low-luminosity quasars was 
addressed by Eisenstein & Loeb (1995), and Loeb (1997). 

Based on the extrapolated quasar luminosity function, we quantify the effects of the first 
generation of low-luminosity quasars. In particular, the UV radiation produced by these quasars 
above the Lyman limit can lead to the reionization of the intergalactic medium (IGM), and to a 
number of observable signatures (Carr, Bond & Arnett 1984). First, the resulting optical depth of 
the universe to electron scattering at z £ 20, r es , damps the microwave background anisotropics 
on scales <; 10° (Efstathiou & Bond 1987; Kamionkowski, Spergel, & Sugiyama 1994; Hu & White 
1996). The recently detected hint for a Doppler peak in ground-based microwave anisotropy 
experiments can already be used to set the constraint r cs 1 (Scott, Silk & White 1995; Bond 
1995). Future satellite experiments (such as MAP or PlanckQ) will be able to probe values of r es 
as small as a few percent if polarization data is gathered (Zaldarriaga, Spergel, & Seljak 1997). 
The quasar UV flux is also absorbed by intergalactic dust, and subsequent thermal emission by 
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this dust at longer wavelengths could introduce a substantial spectral distortion to the cosmic 
microwave background (CMB) radiation. Such a distortion due to population III stars have been 
investigated extensively in the past (Wright 1981; Adams et al. 1989; Bond, Carr & Hogan 
1991; Wright et al. 1994; LH97), and its magnitude is currently constrained by COBE to have 
a Compton y-parameter < 1.5 x 10~ 5 . Finally, direct detection of high-redshift low-luminosity 
quasars will be feasible with the Next Generation Space Telescope, whose proposed sensitivity in 
the 1-3.5/im range is ~lnJy. 

The purpose of this paper is to quantify the above signatures of high-redshift low-luminosity 
quasars, and to compare them with those of early stars. The paper is organized as follows. In § 2, 
we describe the procedure used to find the universal quasar light-curve, and the resulting fits to 
the observed quasar LF at z <J 5. In § 3, we calculate the reionization history of the intergalactic 
medium due to early quasars, and compare it to stellar reionization. In § 4, we discuss the effects 
of the early quasars on the CMB, including the smoothing of its anisotropies and the distortion of 
its spectrum. In § 5, we predict the number counts of high-redshift quasars down to the sensitivity 
limit of the Next Generation Space Telescope. Finally, § 6 summarizes the main conclusions and 
implications of this work. 

2. Quasar Light Curve 

Our goal is to calculate the quasar luminosity function at faint magnitudes, 
(log(LB/LB,©) ^ 11.5) and at high redshifts {z ^ 4.5), outside the regime of current 
observations. For this purpose, we introduce a new model which extrapolates the observed LF of 
bright quasars at redshifts 2.6 < z < 4.5, based on the Press-Schechter theory and the assumption 
that all quasars admit the same universal light-curve in Eddington units. 

The formation rate of dark-matter halos is given by the Press-Schechter formalism. We 
assume that black hole formation is restricted to halos whose virial temperature exceeds 
T viT ^ 10 4 K, in which atomic and Bremsstrahlung cooling allows the gas to sink to the center 
of the potential well. This virial temperature threshold corresponds to a minimum halo mass 
M ha i £ M min = 1O 8 M [(1 + z)/10]" 3 / 2 . We further assume for simplicity that the ratio of final 
black hole mass to halo mass, e = Mbh/Afhalo is the same for all halos. Estimates of black hole 
masses in nearby galaxies (Kormendy et al. 1997; Magorrian et al. 1997) imply that massive 
black holes weigh roughly a fixed fraction of their bulge mass, over a wide range of bulge masses 
1O 9 M <; Mb u ige ^ 2 x 1O 12 M . Since bulges are the oldest stellar components of galaxies, it is 
reasonable to assume that a fixed fraction of the baryonic content of early galaxies ends up in 
the central black hole. Note, however, that this assumption differs from the model presented by 
Haehnelt & Rees (1993), in which e is assumed to decline exponentially for low-mass galaxies. 

The change in the comoving number density of black holes with masses between and 
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Mbh + dMbh, between redshifts z and z + dz, is governed by the derivative, 



d 2 N hh (M hh ,z) 



dM^hdz 



M bh 



1 d dN pB (M,z) 
e dz dM 





M=t~ 1 M h 



(for M bh > eM min ) 
(for M bh < eM min ), 



(1) 



where dN ps /dM is the Press-Schechter mass function. The actual halo formation rate is larger 
than the derivative -^(dN ps /dM), since this derivative includes a negative contribution from 
merging halos. However, at high redshifts collapsed objects are rare, and the merger probability is 
low. We have compared the above expression with the more accurate result for the halo formation 
rate given by Sasaki (1994), and confirmed that the difference in the rates is negligible for the high 
redshifts and halo masses under consideration here. We therefore use equation (||) to describe the 
black hole formation rate at high redshifts. 

We define 4>(L, z)dL to be the number of quasars per unit comoving volume at redshift z, 
with intrinsic luminosities between L and L + dL. The luminosity function, <p(L,z), is linked to 
the black hole formation rate through the light-curves of the quasars. For simplicity, we make 
the minimal assumption that all quasars admit the same light-curve in Eddington units, i.e. their 
luminosity is proportional to their final mass (see §1). Under this assumption, the luminosity L(t) 
of a quasar with a black hole mass Mbh at a time t after its birth can be written as 



L(t) = M hh f(t) = eM halo /(t), (M halo > M min ) 

where f(t) is a function of time. This leads to the expression for the LF, 

, d 2 N hh 
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where t z ^i is the time elapsed between the redshifts z and z' . By integrating over z' , we obtain 

dM bh d 2 N hh 
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or, equivalently, by integrating over Mbh we get 

/■*(*) dt dz 

(j>(L,z) 



f(t) dt 



d 2 N hh 



dMhhdz 



(5) 



M h 



'f(t) 



where / = df '/dt, t* is defined by the condition Mhhf{t+) = L, and z+ is defined through 
the relation = t{z) — t(z+). We found it useful to utilize both of the above expressions for 
(j>(L, z), since a comparison between them can be used to check for the numerical accuracy of the 
calculation. The only unknown quantities are / in equation (Q), and / in equation (||). Therefore, 
either of these two equations can be solved to find the f(t) that would fit best the observational 
data. Since the equations are implicit, we have found their solution iteratively. We have optimized 
our initial guess for f(t) as follows. 
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The function d 2 N\^\ 1 /dM\^dz peaks at a mass scale M p k = eM n \(z), where M n \(z) is the 
non-linear mass-scale at the corresponding redshift. Assuming that this function is sharply 
peaked around M p k, i.e. substituting d 2 N^/dM^dz oc 5(M p k[^]) into equation (Q), yields 

4>(L,z) oc oc , (6) 



\dL/dty 



or 

poo 

ioc J <f>(L',z)dL'. (7) 

In the last step, we assumed that dL/dt < 0, i.e. that the quasar's luminosity decreases 
monotonically with time. Equation (|7|) yields the time t it takes the quasar luminosity to drop 
down to L = Afpk/, and provides a "first guess" solution for t(f) and hence for f(t). The physical 
interpretation of the LF in this case is that it merely reflects the time spent by the quasars in each 
luminosity interval. 

To find the light-curve f(t), we need to specify cj)(L,z). Given the observational data, the 
inferred luminosity function depends on the assumed cosmology and quasar spectrum (the latter 
due to ^-corrections). In this paper, we adopt the "concordance model" of cosmological parameters 
highlighted by Ostriker & Steinhardt (1995), i.e. a flat ACDM model with a slightly tilted power 
spectrum (fl m , £7b, h, cr 8h -i , ra)=(0.35, 0.65, 0.04, 0.65, 0.87, 0.96). Convenient expressions in 
this model for the differential volume element, luminosity distance, and time-redshift relation were 
given in terms of elliptic integrals by Eisenstein (1997). An accurate arithmetic fitting formula for 
the growth function that appears in the Press-Schechter formula, is given by Carroll et al. (1992). 
For the average quasar spectrum, we have used the "median" spectrum derived from observed 
spectra of 47 quasars by Elvis et al. (1994). This spectrum is shown in Figure |l[ 

Pei (1995) provides a fitting formula for the B-band (0.44 ^irn rest-wavelength) luminosity 
function of observed quasars, 4>(L^,z), including the redshift evolution both for a flat and an 
open Qtot = 0.2 cosmology. For conversion to Eddington units, we note that a IMq black hole, 
shining at Eddington luminosity with the median flux distribution of Elvis et al. (1994), has a 
B-band luminosity of 5.7 x 10 3 times the solar B-band luminosity. Pei (1995) assumes power-law 
spectra for the quasars, with a slope a=0.5 in the flat model, and a=1.0 in the open model. Since 
near the B-band the Elvis et al. (1994) template spectrum is very close to an a=l power-law, we 
used the open-universe fitting formulas from Pei (1995), with a rescaling of the luminosities and 
volume elements for our ACDM cosmology. Note that the exact conversion of the intrinsic LF to 
a different cosmology would require a spectrum-dependent ^-correction, as well as rescaling the 
luminosity and the volume element, for the redshift of each individual quasar. However, we have 
found that using the mean redshift in each redshift bin, and assuming the a=1.0 power-law for 
our quasar spectrum, provided the correct LF to within several percent. 

Figure ^ shows the resulting fitting formulas for 4>(Lb, z), together with the observational data 
points for the LF, in the ACDM cosmology at redshifts z=2.6, 3, and 4. Using the fitting formula 
derived by Pei (1995), equation (^) yields an initial guess for f(t). The numerical solution for this 
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guess can be fitted by a parameterized function, and substituted back into equation (|j) or (|5|). 
The parameters of the fit can then be iteratively adjusted so as to optimize the match between 
the resulting 4>(L, z) and the observed quasar LF at z = 2.6, 3, and 4. Using this procedure, we 
have found that an excellent fit is achieved through a simple single-parameter function for the 
light-curve, 

where -£>Edd = 1-4 X 10 38 erg s _1 (Mbh/Af©) is the Eddington luminosity for the finalQ black hole 
mass Mbh- Our model, therefore, has only two free parameters, to and e. The best fit values of 
these parameters are to = 10 5 ' 82 yr and e = 10 -3 ' 2 . The best-fit light-curve is shown in Figure [||, 
while the resulting luminosity functions at z = 2.6, 3, and 4, are shown by the solid lines in 
Figure |2[ 

Our procedure results in fits to the LF comparable in quality to the original ad-hoc parametric 
fits given by Pei (1995). It is rather remarkable that, in the range of observed luminosities, the 
two curves are almost indistinguishable, even though at fainter magnitudes our curves deviate 
from the Pei (1995) fits, and predict many more quasars than the extrapolation of those fitting 
formulae. Another gratifying feature of our approach is that - with the underlying assumption 
that the flat portion of the light -curve corresponds to the Eddington luminosity limit of the final 
black hole - it results in a prediction for the black hole to halo mass, M^/M^io = 10~ 3 ' 2 , or black 
hole to gas mass Mbh/M gas = 10 _3 ' 2 $7 m /T2b = 5.4 x 10 -3 . This is close to the value derived from 
observational data, 6 x 1CP 3 , for the bulges of nearby galaxies (Kormendy et al. 1997; Magorrian 
et al. 1997), which might have formed at about the same time as quasars. There are, however, 
two shortcomings of our model. First, since it is based on the Press-Schechter formalism, it 
predicts that the number of bright quasars (those with masses above M p k, or Lb ^ 10 13 Lb,©) 
is steadily increasing down to z = 0. Hence, by construction, it cannot explain the decrease in 
the number density of bright quasars at low redshifts (z ^ 2). This decrease might be associated 
with the decline in merger rates at low redshifts (Carlberg 1990); or it might also be a result of 
the expulsion of cold gas from galaxies through supernovae-driven winds due to the peak in the 
star formation activity around z ~ 2 (Madau 1997). However, since our predictions focus on 
high redshifts (z ^ 5), we ignore this shortcoming of the model. Second, the best-fit light-curve 
predicts a quasar lifetime of only ~ 6.6 x 10 5 years, about two orders of magnitude shorter than 
the Eddington time at the usually assumed value of the radiative efficiency, 10%. The average 
radiative efficiency over the growth history of the black hole, is ~0.1%, two orders of magnitude 
smaller than the value for thin disks (Frank et al. 1992). This may imply that quasar black 
holes accumulate most of their mass in an accretion flow with a low radiative efficiency, such as 



4 During the initial growth of the black hole, its mass and hence its Eddington limit, are smaller. This is not in 
conflict with our prescribed light -curve, as long as most of the light is emitted during a late phase of duration to, 
during which the black hole shines close to the Eddington limit of its final mass, Afbh- The inferred low value of 
to implies that much of the black hole mass is accumulated via accretion with a low radiative efficiency (e.g., an 
advection-dominated accretion flow, cf. Narayan 1996 and references therein). 
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advection dominated flows (Narayan 1996, and references therein). Although to is an unusually 
short lifetime for the bright phase of quasars, it is not in conflict with existing observations. 
Furthermore, increasing the lifetime of quasars would enhance all of the effects that we describe 
below. In this sense, we are being conservative by adopting our best-fit light-curve. 

The evolution predicted by equation (||) for the quasar LF between redshifts 2.6 < z < 18, is 
shown in Figure |4|. The break at Lb ~ 1O 1o Lb,0 is introduced by our imposed low-mass cutoff 
-^haio > 1O 8 M [(1 + z)/10]~ 3 / 2 . In Figure [f| we show the evolution of the number density of 
quasars brighter than absolute blue mag nitude M B = -27.5 and -22.5 (L B = 10 13 - 2 L B ,© and 
1O 11 ' 2 Lb,0, respectively). The solid lines in this figure show the number densities predicted by 
our model, while the dotted lines show the corresponding number densities based on Pei's (1995) 
fitting formulas. This figure demonstrates that although our model agrees with the observational 
data between 2.6 < z < 4.5, it predicts a substantially larger abundance of quasars at higher 
redshifts and fainter magnitudes, than the formal extrapolations of the fitting formulae given by 
Pei (1995). 

3. Reionization 

In this section, we derive the reionization history of the IGM in the presence of early quasars. 
The method of calculation follows closely the one introduced in HL97 for the case of stellar 
reionization. We refer the reader to this reference for more details. 

Each ionizing source creates an expanding Strbmgren sphere around itself. Reionization 
is complete when the separate HII regions overlap (Arons & Wingert 1972). Since the neutral 
fraction inside each HII region is small (~ 10 -6 ) and the mass- weighted dumpiness of the IGM 
is still small at z ;> 10 (Gnedin & Ostriker 1997), the average ionized fraction at any redshift 
is simply given by the filling factor of the HII regions. The reionization history depends on the 
formation rate of quasars, the light-curve of ionizing photons from quasars, and the recombination 
rate in the IGM, which are all functions of redshift. 

The production rate of ionizing photons per quasar follows from the quasar spectrum and 
light-curve shown in Figures [l] and ^. The production rate per solar mass of the black hole, 
dNE>i3.QeV /dt = 6.6 x 10 47 Mbh f{t) photons Mq 1 s _1 , is shown as the solid line in Figure ||. For 
comparison, we also show the analogous rate per stellar mass from the composite stellar spectrum 
in our model of a zero-metallicity starburst. Note that although the quasar rate is initially higher 
than the stellar rate by two orders of magnitude, it drops much more rapidly at later times. The 
total number of ionizing photons produced by a quasar is ~ 1.3 x 10 61 M bh /M , only an order of 
magnitude higher than that of stars, ~ 1.3 x 10 60 M stars /M0. 

We assume that recombinations take place only in the homogeneous IGM. This includes the 
implicit assumption that all of the ionizing photons escape from the neighborhood of the quasar. 
While this would be a poor assumption for stellar radiation, it is justified for the much more 
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intense (and harder) quasar emission. A simple estimate can be used to show that the quasar host 
would be ionized shortly after the quasar turns on. Indeed, the known quasars do not show any 
sign of associated (galactic) HI absorption in their spectrum, while stellar spectra of high redshift 
galaxies often do. An additional assumption in obtaining the recombination rate is that the IGM 
is homogeneous, i.e. the clumping factor for ionized hydrogen is unity, Chii = (^hii) / ( n Hii) 2 = 1- 
This is justified in numerical simulations, which show that Chii rises significantly above unity only 
when the collapsed fraction of baryons is large, i.e. at z <^ 10 (see Fig. 2 in Gnedin & Ostriker 
1997) 

The evolution of a cosmological ionization front for a time-dependent source is governed by 
the equation, 

"» - H{t)r ) = it? (^ (t » - r*-"*) • (9) 

where «h is the neutral hydrogen number density in the IGM, r\ is the physical (i.e. not comoving) 
radius of the ionization front, H(t) is the Hubble expansion rate at cosmic time t, dN^/dt is the 
ionizing photon rate, and «b = 2.6 x 10~ 13 cm 3 s _1 is the recombination coefficient of neutral 
hydrogen at T = 10 4 K (see HL97, or Shapiro & Giroux 1987 for more details). The solution of 
this equation yields the ionized volume per source, 4/3nrf. A summation over all sources provides 
the total ionized volume within the IGM. 

The results for the stellar case depend on the assumed star formation efficiency, which we 
normalize based on the observed metallicity of the IGM at z ~ 3 (HL97). In individual Lya 
clouds with JVhi *5 10 14 cm -2 , where C and/or Si are detected, the inferred metallicities (in solar 
units) are ~ lCT 2 5 for C, and ~ 1(T 2 for Si (Songaila & Cowie 1996, Tytler et al. 1995). There 
appears to be no trend of a decreasing mean C/H or Si/H with decreasing column density, across 
six orders of magnitude in JVhi. Hence, a plausible assumption is that the metals are universally 
mixed into the IGM prior to z = 3, so that the metallicities of the low column density absorbers 
are the same as that found in systems with iVm *5 10 14 cm -2 . However, as pointed out by 
Songaila (1997), it is also possible that there are no metals in the low column density systems; in 
that case, a conservative lower limit for the average metallicity is Zjqm = 1O _3 Z0. The resulting 
range of IGM metallicities, 1CT 3 < (^igm/Z©) < 1CP 2 , translates to the corresponding range 
0.017 < / s tar = Mstax/Mgag < 0.17 for the star formation efficiency. 

Figure ^ summarizes the resulting reionization histories of stars or quasars in our ACDM 
cosmology. The results for stars are shown in two cases, one with Ziqm = 10 _2 Zq and 
the other with Ziqm = 1O _3 Z0, to bracket the metallicity range mentioned above. The 
assumed IGM metallicity does not affect the quasar calculation. The upper left panel of 
Figure [7| shows the fraction of baryons which reside in collapsed objects with a mass above 
M min = 1O 8 M [(1 + z)/10]- 3 / 2 . The upper right panel shows the evolution of the average flux, 
J21, at the local Lyman limit frequency, in units of 10 -21 erg s _1 cm -2 Hz -1 sr~ 4 for quasars 
(solid line) and stars (dashed line for Ziqm = 1O _2 Z ; and dotted line for Ziqm = 1O _3 Z0). The 
equations used to obtain J21 are the same as described in LH97, except for the changes in the line 
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element, time-redshift relation, and baryonic collapsed fraction, caused by our switch to a ACDM 
cosmology. The stellar mass function is taken to have the Scalo (1986) form. We have assumed 
that the flux between the local Lya and Lyman-limit frequencies vanishes before reionization, due 
to the large Gunn-Peterson optical depth. Although this is a good approximation for the present 
purpose, we note that the actual spectrum in this frequency range is expected to have a sawtooth 
shape, as described in Haiman, Rees & Loeb (1997). After reionization J21 rises rapidly, and by 
a redshift z = 3 the stellar J21 reaches a value in the range 0.04-0.4, while quasars contribute 
J 2 i ~ 0.1. These values are consistent with the background flux level inferred from the proximity 
effect around quasars (Bechtold 1994). 

The lower left panel of Figure 7 shows the resulting evolution of the ionized fraction of 
hydrogen, -Fhii- The dashed curves indicate that for the high end of the allowed metallicity range, 
stars ionize the IGM by a redshift z ~ 13; while the dotted curve shows that for the lowest 
metallicity, reionization is delayed until z = 9. (Note that for Z\qm = O.OIZq, HL97 obtained 
stellar reionization at a redshift z = 18 in standard CDM. The change to z = 13 here results 
from our switch to a ACDM cosmology.) The solid curve shows that quasars reionize the IGM 
between these two redshifts, at z = 11.5. This result can be understood simply in terms of the 
total number of ionizing photons produced in each case per unit halo mass: the relative ratios of 
this number in the three cases are 1 -j- 0.37 -^0.1, respectively. 

Although the reionization redshift can be inferred from the corresponding damping of CMB 
anisotropies (see below), it is not possible to determine whether reionization was caused by stars 
or quasars, based on the measurement of this damping alone. However, the distinction between 
the two possibilities could be made if the reionization redshift of Hell is also measured. If a 
source at a redshift prior to reionization is discovered, the location of the Gunn-Peterson trough 
in its spectrum would yield the reionization redshift, while the shape of the damping wing of the 
trough can be used to infer the Gunn-Peterson optical depth (Miralda-Escude 1997). Since the 
ionization threshold of Hell is four times higher than that of HI, the ionization of Hell requires 
a much harder spectrum. The relevant measure is rj = (dNE>54AeV / dt) / (dNE>i3.6eV /dt), the 
ratio of photon fluxes above 54.4eV and 13.6eV. Equation (||) reveals that if recombinations are 
negligible, and if 77 < nn e /nn ~ 0.08, then the Hell/Helll ionization front will lag behind the 
HI/HII-ionization front (see Miralda-Escude & Rees 1994, for a discussion). If recombinations 
are important (redshifts z 15), then the lag is increased further, since Hell recombines 5.5 
times faster than HI. In our models, the Hell reionization history can be evaluated by replacing 
nu with rc-He, and dNE>iz.&e\ / dt with dNE>54Aev/dt, and multiplying «b by a factor of 5.5, in 
equation (^). 

We find that for our quasar spectrum, HI and Hell reionization occurs almost simultaneously 
(the Hell reionization redshift is z = 11.8), but for our stellar spectrum, Hell reionization is 
never achieved. These results are explained by fact that for the quasar spectral template of 
Elvis et al. (1994), n = 0.09; while for the stellar spectral template, r\ ~ 10~ 6 , i.e. there are almost 
no Hell-ionizing photons, because most of these photons are absorbed already in the stellar 
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atmospheres (see Fig. 3 in HL97). It therefore seems likely that intergalactic Hell was reionized 
by quasars rather than by stars. Thus, the reionization history of Hell provides an ideal probe for 
the formation history of the first quasars. 

Our results indicate that the reionization redshifts of HI and Hell should not be very 
different, unless we have substantially overestimated the emission from high-redshift quasars, or 
underestimated the abundance of high mass stars. Reimers et al. (1997) have recently claimed 
a detection of the Hell reionization epoch at z = 3. Their claim is based on the lack of any 
detected flux above the redshifted Hell Lya frequency, although the flux at the corresponding 
redshifted HI Lya frequencies is transmitted without any absorption. This observation, however, 
only places the limit THe ^ few on the Hell Gunn-Peterson optical depth, rather than the th c <; 
few thousand, required to prove Hell reionization. In fact, Miralda-Escude (1997) argued that 
these observational results are likely to have been caused by spatial variations in the gas density, 
and fluctuations of the background radiation intensity above the Hell edge, rather than the fact 
that the Helll bubbles have not yet overlapped. In summary, the interpretation of this, and other 
recent detections of Hell absorption (Jacobsen et al. 1994; Tytler et al. 1995; Davidsen et al. 1996; 
Hogan et al 1997) are still controversial, in that these measurements are consistent with the 
absorption arising either from non-overlapping Helll bubbles, or from other effects associated with 
intervening discrete systems (see, e.g. Madau &: Meiksin 1994; Giroux & Shull 1997). 

4. Signatures Imprinted on the Cosmic Microwave Background 

Reionization results in a reduction of the temperature anisotropies of the cosmic microwave 
background (CMB). The free electrons released during this epoch increase the optical depth 
of the universe to electron scattering (r es ), and damp the CMB anisotropies on scales below 
the size of the horizon at that time, ^ 10° (Efstathiou & Bond 1987; Kamionkowski, Spergel, 
h Sugiyama 1994). The electron scattering optical depth in our models is approximately 
T es (z) w 0.053f2b^ / dz\J\ + zFun(z), where Fmi(z) is the filling factor of the HII regions. Hu 
& White (1997) provide a fitting formula for the damping factor of the CMB power-spectrum, 
Rg, as a function of the index I in the spherical harmonic decomposition of the microwave sky 
[the angular scale corresponding to a given ^-rnode is ~ 1° x (^/200) -1 ]. The damping factor 
is uniquely related to the evolution of the optical depth for electron scattering as a function of 
redshift. 

The lower right panel of Figure shows both T es (z) and R\ for the stellar and quasar 
reionization histories. The electron scattering optical depth is 3-5% for the stellar case, and 
~4% for the quasar case. The corresponding damping factors for the power spectrum of CMB 
anisotropies are 6-10%, and ~8%, respectively. Although the amplitude of these damping factors 
is small, they are within the proposed sensitivities of the future MAP and Planck satellite 
experiments, if data on both temperature and polarization anisotropies of the CMB will be 
gathered (see Table 2 in Zaldarriaga et al. 1997). 
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The UV emission from early quasars could also distort the spectrum of the CMB, after 
being reprocessed by dust. An early epoch of star formation and metal enrichment is inevitably 
accompanied by the formation of dust in supernova shells. This dust would absorb the UV flux 
from both stars and quasars, and re-emit it at longer wavelengths. The re-emitted radiation is 
added to the CMB spectrum, and introduces a deviation from its perfect blackbody shape. A 
distortion of this type due to population III stars alone has been considered by several authors 
(Wright et al. 1994; Bond, Carr k Hogan 1991; Adams et al. 1989; Wright 1981). 

In a previous paper (LH97), we have calculated the distortion amplitude from early stars 
in a standard CDM cosmology, using the Press-Schechter formalism and normalizing the star 
formation efficiency in collapsed objects based on the observed metallicity of the IGM. The 
essential assumptions in this calculation were: (1) the initial mass function (IMF) of stars at 
high-redshifts is the same as observed locally (Scalo 1986), and (2) each type II supernova yields 
0.3Mq of dust, which gets uniformly distributed within the intergalactic medium. We further 
assumed that the dust is in thermal equilibrium with the total (CMB + stellar) radiation field, 
that its absorption follows the wavelength-dependent opacity of Galactic dust (Mathis 1990), and 
that it radiates as a blackbody at its equilibrium temperature. We have found that allowing two 
dust components made of graphite and silicate (Draine & Lee 1984) with independent equilibrium 
temperatures, does not change our results by more than a few percent. For a range of possible 
parameter values within standard CDM, we found that the opacity of intergalactic dust to infrared 
sources at redshifts z ^ 10 is significant, Tdust = (0.1-1), and that this dust distorts the microwave 
background spectrum by a Compton y-parameter in the range (0.06-6) x 10~ 5 . 

Here, we repeat the LH97 calculations in a ACDM cosmology, and extend them by adding the 
UV flux from quasars to the total radiation field. Our new results are presented in Figure || The 
top panel of this figure shows the resulting total spectrum of the radiation background (CMB + 
direct stellar and quasar emission + dust emission) at z = 3. The spectrum describes a minimum 
level of CMB distortion down to z = 3; more distortion could be added between < z < 3 by dust 
and radiation from galaxies. However, the additional emission at < z < 3 does not relate directly 
to the reionization epoch and is therefore ignored here. The dashed lines show the spectrum 
for the stellar case, for Z\qm = 1O~ 2 Z (top curve) and Zjqm = 1O~ 3 Z (bottom curve). The 
pronounced peak at a wavelength of ~ 1/im is the sum of the direct starlight from high-z stars, 
redshifted to z = 3. The deviation from the pure Tcmb = 2.728(1 + z)K blackbody shape (shown 
by the dotted curve) are visible in this figure below a wavelength of 0.02cm. The corresponding 
distortions in the COBE range is quantified by the Compton y-parameter, defined as 



1 



IduN u 
I duN , v 



(10) 



where the integrals are evaluated over the FIRAS frequency range of 60-600 GHz (Fixsen et al. 
1996), N u is the comoving number density of photons with a comoving frequency v in the total 
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radiation field, and Nq >v reflects the unperturbed CMB, 

The bottom panel of Figure || shows the redshift evolution of the y-parameter, as well as the dust 
temperature, T^ust- For the stellar cases considered above, we obtain 1.3 x 1CT 7 < y c < 1.2 x 1CT 5 
at z = 3 (dashed lines). Since both the dust opacity and the stellar radiation background 
are independently proportional to the star-formation efficiency /star, the spectral distortion is 
proportional to / s 2 tar (see LH97). For the high end of the allowed metallicity range, the distortion 
at z = 3 is just below the upper limit set by COBE, y = 1.5 x 10 -5 (Fixsen et al. 1996); for lower 
metallicities, the y-parameter drops rapidly below this value. The dust temperature remains in 
all cases close the the CMB temperature (shown by the dotted line). 

The solid lines in Figure || show how these results change when the radiation from the early 
quasars is added. In particular, the top panel shows that the broad emission spectrum of quasars 
(cf. Fig. 1) adds a background radiation field over a much broader range of wavelengths than 
stellar emission. But since dust is effective in absorbing the quasar flux only in a relatively narrow 
range around ~ l^m, the distortion in the COBE regime is not significantly altered. In particular, 
we find 4.1 x 1(T 6 < y c < 2.0 x 10~ 5 at z = 3 for the range 1O" 3 Z < Z IGM < 1(T 3 Z . It is 
interesting to note that the long- wavelength emission from high redshift quasars in the COBE 
wavelength range, results in a minimum distortion of y c ~ 3.4 x 10~ 6 . Depending on Z\qm, this 
amounts to a fraction ~ 17-82% of the total y-parameter values quoted before. This minimal 
level of distortion results simply from the cumulative far-infrared flux of early quasars, and should 
be present even in the absence of any intergalactic dust. It can be avoided only if early quasars 
have much less emission longward of ~ 100/um, relative to the spectrum shown in Figure 1. 



5. Number Counts of Faint Quasars 

Next we examine the feasibility of direct detection of the early population of low-luminosity 
quasars. Although any such detection is beyond the sensitivities of current observing programs, 
future instruments such as the Space Infrared Telescope Facility (SIRTF) or the Next Generation 
Space Telescope (NGST), will achieve the required sensitivity. NGST, which is scheduled for 
launch a decade from now, is expected to reach a sensitivity of ~ 1 nJy in the wavelength range 
1-3. 5/um (Mather Sz Stockman 1996; see also htt p : / / ngst . gsf c . nasa. gov] ) . The number of sources 



that NGST would detect at redshifts z > z m i n , per unit solid angle, with observed flux between 
F v and F v + dF u (averaged over the observed wavelength range 1-3.5/xm), is given by 



(F u ,z min )= [°° dz(^-)n c (z,F u ), (12) 



dN 

d£ldF v x Jz min \dzdflj 

where dV c /dzd£l is the comoving volume element per unit redshift per unit solid angle, and 
n c (z,F u ) is the comoving number density of objects at redshift z, whose flux at z = is observed 
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to be between F v and F u + dF v . Based on the halo formation rate d 2 N ps /dMdz [cf. Eq. (|])], 
given by a sum over halos of different ages that exist at each redshift, 



n c (z,F u 



00 HM d 2 N 



z 



dF v 



dMdz 



n r is 



(13) 



M{z,z',F v ),z' 



where the factor dM/dF u converts the number density per unit mass interval to number density 
per unit flux interval, based on equation (]l~4|). The assumptions implicit in this expression are that 
the simple derivative d 2 N ps /dMdz gives the halo formation rate [see discussion following eq. (EI)], 
and that the halos do not grow in mass during the active phase of the emission sources they host 
(a good assumption for the short-lived quasar activity) . Relaxing any of these assumptions would 
increase our predicted number counts. Note that, in HL97 we had effectively identified z = z' in 
the kernel of equation (|i~3|), when estimating the number counts of star clusters. This might be 
justified only for stellar sources that live for a Hubble time or longer. However, in this case there 
is a danger of double counting merging halos in neighboring redshift bins (although at sufficiently 
high redshifts objects are exponentially rare and their merging probability is small). We intend to 
address this issue more rigorously using the excursion-set formalism (Bond et al. 1991; Lacey & 
Cole 1993, 1994), in a future publication. 

The average flux F u (z, M) from a halo at redshift z with mass M, that was formed at z' > z 
equals 

10 32 eM 



F u (z,z') 



dvj{v, t z 



nJy. 
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d L (z) 2 Az, 

Here the halo mass M is in units of M Q ; z^ m j n and ^ max are the redshifted frequencies (in units of 
Hz) corresponding to the observed wavelengths of 3.5^m and 1/xm, respectively; Avq = 2.14 x 10 14 
Hz is the width of the NGST frequency band; d^{z) is the luminosity distance in cm; e = 10 -32 
for quasars or e = 0.17fib/^m = 0.02 for stars, is the fraction of halo mass in the form of the 
central black hole or stars, respectively; j{v,t z ^ z i) is the template luminosity per unit (stellar or 
black hole) mass for either stars or quasars, in units of erg sec -1 Hz -1 sr _1 Mq 1 ; and t z ^ z > is the 
time elapsed between redshifts z' and z. The time dependence of j(y) is obtained from our fitted 
light-curve function f(t) in equation (g) for quasars, and from tabulated evolutionary tracks for 
stars (Schaller et al. 1992; see HL97 for details). We suppressed the emergent flux from all objects 
at emission frequencies above Lya before reionization occurs (assuming sudden reionization at 
the redshift z = 13 or 9 for stars, and z = 11.5 for quasars), since the Gunn-Peterson optical 
depth of HI at these frequencies is exceedingly high prior to reionization. 

Figure |9| shows the predicted number counts normalized to the field of view of SIRTF (which 
is scheduled for launch long before NGST), 5' x 5'; the planned field-of-view of NGST is similar, 
4' x 4'. Figure ^ shows separately the number per logarithmic flux interval of all objects with 
z > 5 (thin lines), and with z > 10 (thick lines). In general, the number of detectable sources is 
high. With its expected sensitivity, NGST will be able to probe about 10 3 quasars at z > 10, 
and 10 35 quasars at z > 5 (solid curves) per field of view. The bright-end tail of the number 
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counts approximately follows the power law dN/dF u oc F~ 2 ^. This slope is too shallow for 
foreground galaxies to significantly increase the number counts due to the amplification bias from 
gravitational lensing. The top solid curve also shows that the number of faint quasars starts to 
turn over below ~ 1 nJy, as a result of the cutoff in the minimum halo mass that allows atomic line 
cooling. The average angular separation between the quasars at z > 10 would be ~10", well above 
the angular resolution limit of 0.06" planned for NGST. For comparison, we show in Figure ^ 
the corresponding number counts of "star-clusters", i.e. assuming that each halo shines due to a 
starburst that converts a fraction 0.017-0.17 of the gas into stars. The dashed lines indicate that 
NGST might detect 10 L5 -10 3 - 5 star-clusters at z > 10 per field of view, and 10 3 - 5 -10 4 5 clusters 
at z > 5. Unlike quasars, star clusters could in principle be resolved, if they extend over a scale 
comparable to the virial radius of the collapsed halo (Haiman & Loeb 1997b). 

6. Conclusions 

We extrapolated the quasar luminosity function to faint magnitudes and high redshifts, based 
on the assumptions that the halo formation rate follows the Press-Schechter theory, and that 
the quasar light-curve scales linearly with the halo mass, and is otherwise a universal function 
of time. We demonstrated that a universal light-curve of the form, L(t) = L^dd exp(— t/to), 
with to = 6.6 x 10 5 yr and i^Edd being the Eddington luminosity, provides an excellent fit to 
observational data for the quasar LF between 2.6 < z < 4.5 provided that the final black hole mass 
is a fraction ~ 5.4 x 10~ 3 of the baryonic mass in each object. This fraction is surprisingly close 
to the roughly universal ratio between black hole and bulge masses observed in nearby galaxies 
(Magorrian et al. 1997). Our extrapolation of the quasar LF to high redshifts and low luminosities 
as presented in Figure |||, can be viewed as a natural by-product of this physically motivated fit. 

Given the Press-Schechter prediction for the formation rate of halos, and the fact that to is 
much shorter than the Hubble time, our prediction for the reionization history depends only on the 
integral of the inferred light-curve, namely on the mean radiative efficiency of the material that 
made the black hole. The low value inferred for to implies an average radiative efficiency of only 
~ 0.1%. This rather low estimate for the efficiency is smaller by two orders of magnitude than the 
value expected in thin disk accretion; any increase in its value would only enhance the effects we 
predict. Despite our conservative approach, we predict an extensive population of faint quasars at 
high redshifts. For illustration, we show in Figure ||| the expected evolution in the number density 
of bright quasars with absolute B magnitude Mb < —27.5, along with that of quasars which are 
two orders of magnitude fainter. Although the density of bright quasars in our models declines 
rapidly beyond redshift z ~ 3 in good agreement with existing observations at 2.6 < z < 4.5, we 
predict an abundant population of faint quasars which extends out to redshifts z p> 10, and has 
not yet been directly observed. 

We find that, in general, the effects of early quasars are comparable to those of early stars. 
However, stellar effects scale with the average IGM metallicity that we used to normalize the star 
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formation efficiency. If the observed values of the C /H ratio in the individual Lya forest clouds 
are representative of the universal average C/H ratio, than this ratio implies an IGM metallicity 
of ~ 0.5-1% solar. However, if the overall carbon abundance is patchy, because it is not well 
mixed with the rest of the intergalactic gas (as argued by Gnedin 1997), then the true average 
metallicity may be an order of magnitude lower, 0.1% solar (Songaila 1997). We have found that 
at the high end of the allowed metallicity range, the effects of stars are stronger than that of 
quasars, while at the low end of this range, quasars are more effective. In particular, quasars 
reionize the universe at z ~ 11.5, while stars at 9 ^ z ^ 13. The Compton y-parameter for the 
spectral distortion of the CMB is expected to be 1.3 x 10 -7 < y c < 1.2 x 10~ 5 from starlight 
alone, and 4.1 x 10 -6 < y c < 2.0 x 10~ 5 if the flux of quasars is added. The direct far-infrared 
flux from quasars yields a lower limit to the distortion, y c ^ 3.4 x 10~ 6 , even in the absence of any 
intergalactic dust. The expected number of z > 10 quasars per square arcminute, brighter than 
InJy in the wavelength range 1-3.5/um, is ~ 40, while the corresponding number for star clusters 
is ~ 1-100. Both of these numbers are large enough to show as a distinct population on a typical 
image by the forthcoming Next Generation Space Telescope. 

We thank J. McDowell, Y. Pei, and W. Zheng for data on quasar spectra and the quasar 
luminosity function, and L. Hui, J. Miralda-Escude, R. Narayan, M. Rees, E. Turner, and A. 
Songaila for useful discussions. This work was supported in part by the NASA ATP grant 
NAG5-3085, and the Harvard Milton fund. 
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Fig. 1. — Comparison of the observed spectral fluxes from a typical quasar and a starburst in a 
source located at redshift z s = 10. The star-formation efficiency, / s tar = 0.17, was fixed so as to 
produce an average of 1% solar metallicity in the IGM at z = 3. The black hole formation efficiency, 
Mbh/Mhaio = 10~ 3 ' 2 , was found by our fitting procedure to the observed quasar luminosity function 
(see text for details). 
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Fig. 2. — The observed and extrapolated evolution of the quasar luminosity function in a ACDM 
cosmology. The data is taken from Pei (1995), and the dashed lines show the parametric fitting 
function from this reference. The solid lines show our fits based on the Press-Schechter formalism 
and the universal light-curve shown in Fig. 3. 
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Fig. 3. — The best-fit quasar light-curve in Eddington units, that results in the luminosity 
functions shown by the solid lines in Figure |2[ The Eddington time was not used to rescale the 
time axis, because it does not depend on the black hole mass, and is a constant for a universal 
radiative efficiency. 
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Fig. 4. — The extrapolated quasar luminosity function at high redshifts, using the light -curve of 
Figure |. The break at L B ~ 10 10 L B . is a result of the low mass cutoff for the dark matter halo, 
M m i n = 1O 8 M0[(1 + z)/1O] _3//2 , that defines the threshold for efficient (atomic) cooling of virialized 
gas. 
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Fig. 5. — The redshift evolution of the number density of quasars. Solid lines show the predictions 
from our model, while the dashed lines show the formal extrapolation of the fitting formulae 
obtained by Pei (1995). 




Fig. 6. — Evolution of the production rate of ionizing photons, per unit black hole mass in quasars 
(solid curve), and per unit stellar mass in a starburst (dashed curve). 
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Fig. 7. — Reionization history from quasars (solid lines) and early stars (dashed [Zjqm = 1O _2 Z0] 
and dotted [Ziqm = 1O _3 Z0] lines). The different panels show (i) the collapsed fraction 
of baryons; (ii) the evolution of the flux at the local Lyman limit frequency, in units of 
1CT 21 erg s _1 cm~ 2 Hz -1 sr _1 ; (hi) the volume filling factor of ionized regions; and (iv) the optical 
depth to electron scattering, r es , and the corresponding damping factor, R|, for the power-spectrum 
decomposition of CMB anisotropies as a function of the spherical harmonic index i. 
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Fig. 8. — Upper panel: The full (CMB + direct light + dust emission) background flux at redshift 
z = 3 in the presence of early stars alone (dashed lines) and in the presence of both stars and 
quasars (solid lines). The pairs of dashed and solid lines bracket the possible range of the IGM 
metallicity, 10 _3 Z Q < Zjqm < 1O~ 2 Z . Lower panel: The corresponding redshift evolution of the 
dust temperature, and the Compton y-parameter. The dotted curve shows the CMB temperature, 
Tcmb = 2.728(1 + z)K. 



- 28 - 



|\ wavelength range: 1 3.5 /xm 




12 3 4 

log (FynJy) 



Fig. 9. — The number per logarithmic flux interval of high-redshift objects which could be probed 
by future space telescopes, in the wavelength range of 1-3.5/im. We assume sudden reionization at 
z=11.5 for quasars, and at z=9 or 13 for stars. The vertical dashed lines show the expected 
sensitivities of the Space Infrared Telescope Facility (SIRTF) and the Next Generation Space 
Telescope (NGST). The thick lines, labeled "10", correspond to objects located at redshifts z > 10, 
and the thin lines, labeled "5", correspond to objects located at redshifts z > 5. 



